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Abstract 

Three-body systems of scalar bosons are covariantly described in the framework 
of relativistic constraint dynamics. With help of a change of variables followed 
by a change of wave function, two redundant degrees of freedom get eliminated 
and the mass-shell constraints can be reduced to a three-dimensional eigenvalue 
problem. In general, the reduced equation obtained by this procedure involves 
the spectral parameter in a nonconventional manner, but for three equal masses 
a drastic simplification arises at the first post-Galilean order: the reduced wave 
equation becomes a conventional eigenvalue problem that we treat perturbatively, 
computing a first order correction beyond the nonrelativistic limit. The harmonic 
interaction is displayed toy model. 



1 Introduction, Basic Equations 

A relativistic system of mutually interacting particles can be described, in a man- 
ifestly covariant way, by mass-shell constraints. These constraints determine the 
evolution of a wave function which depends on four-dimensional arguments [1-3]. 
The price paid for covariance is the presence of redundant degrees of freedom, just 
like in the Bethe-Salpeter (BS) approach. 

For two-body systems, the extra degree of freedom is trivially factorized out. More- 
over, in that case, the contact with the BS equation was established [4] . 
In contrast, for three or more particles it is difficult to find an interaction term such 
that the mass-shell constraints are compatible among themselves, respect Poincare 
invariance, reproduces free-body motion when this term vanishes, and allow for 
eliminating the redundant degrees of freedom. 

In this paper we focus on the case of three spinless particles. We thus consider 
three Klein-Gordon equations coupled by a mutual interaction which should be 
either derived from the underlaying field theory (QCD for instance) or motivated 
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by phenomenological considerations. Our basic equations are 

2K a $ = (p 2 a + 2W) $ = m 2 a c 2 $, a, b =1,2, 3 (1) 

for a wave function with three four-dimensional arguments (either q a or p b accord- 
ing to the representation used). The relativistic "potential" W must be Poincare 
invariant and chosen such that the equations above are mutually compatible. 
Remark 

In principle it seems that one could also consider more general equations involving 
three distinct relativistic "potentials" W\, W2, W3. But, even if we leave aside the 
problem of mutual compatibility which would become more complicated, in this 
more general case there is no evidence that the superfluous degrees of freedom 
could be eliminated at all. It is therefore natural to focus on the simple class of 
models characterized by W a = W . This choice is reminiscent of what is currently 
done in the two-body case, where using the same interaction function in both wave 
equations is general enough to accommodate most realistic situations [5] [6]. 
For three particles, the assumption of a single interaction function in eqs (1) will be 
justified a posteriori by its efficiency in the task of reducing the degrees of freedom. 
Of course in a future work, a justification by a contact with field theory would be 
desirable. 

It should be emphasized that having W a = W by no means forbids to take into 
account differences in the couplings that concern each particle. In most systems of 
practical interest W is a sum of three terms; each one of these three terms, although 
it is not strictly binary (due to the three-particle forces automatically included) still 
carries some two-body input (the three-body forces being of higher order). 

This will be seen for instance in the covariant harmonic potential of Section 
4, where the potential given by equation (51) includes three distinct spring con- 
stants permitting to implement a different interaction law inside each cluster; for 
electromagnetic interactions different charges could be handled in a similar way. 

Poincare algebra is realized in the same manner as for the free-particle case, say 
P = p\ +P2 +P3 and M = J2Qa^Pa- It is convenient to introduce relative variables 
with the "heliocentric" notation: relative-particle indices are A,B = 2,3. We define 
the four-vectors 

P 

z A = qi~ qA, VB = -g- - PB (2) 

Their transverse parts are za,Vb- Tilde denotes the projection orthogonal to P, for 
instance za = za — {za ■ P)P/P 2 , etc. 

With help of the identity 

3 & 2 = P 2 + D + 6P 2 Z (3) 

where 

D = 6(y 2 2 + yl + m- m) (4) 
s = (P 2 r 2 [(y 2 • P? + (y 3 ■ P? + (y2 ■ P)(m ■ P)\ (5) 
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the sum of equations (1) yields a dynamical equation involving the potential. On 
the other hand, the differences of equations (1) take on the purely kinematic form 

(Pi - Pa)(pi + Pa) $ = 2u A c 2 $ (6) 

where the half squared-mass differences are u A = 7^( m i ~~ 171 a) ■ For the sake of 
compatibility we require that W commutes with both products (p\ — Pa)(pi + Pa)- 

In order to achieve the elimination of two degrees of freedom, we have proposed [7] 
a quadratic change of variables in momentum space, say p a ► p'f, or equivalently 

P, y A i-> P, 

This transformation can be characterized as a redefinition of the relative energies 
such that 

(pi + Pa) • (Pi -p A ) = P- (p'i - Pa) (7) 
and by the requirement that it leaves P, yi and ^3 unchanged, that is 

P' = P, y' A = y A (8) 

Clearly equation (7) determines in closed form the longitudinal pieces of p'^Pz (resp. 
y^y's) in terms of all the primitive variables p a (resp. P,yA) [8]. 
Of course, we define the new relative momenta y' A as linearly related with the p'^s 
through a formula similar to (2), namely 

Va = ^~Pa 

Note that our transformation preserves Poincare invariance; as a result the genera- 
tors of spacetime displacements have the usual form also in terms of p' a . 
Naturally, this procedure gives rise to new configuration variables z' A . In general, 
the new variables z' A ,y' B are referred to as reducible. 

It is noteworthy that instead of $, we can equivalently use a new wave function 

' ' D(P'1,P2,P3) 

Accordingly, the operators K and W are mapped to H and V respectively, 

H = \J\ 1,2 K\J\- 1 ' 2 , V = \J\ 1/2 W\J\- 1/2 (10) 

In contrast to Z2,z^, the operators z' A are "formally hermitian" (i.e. symmetric) 
with respect to the Hilbert space L 2 (R 12 , d 12 p'). In other words, setting 
d 12 p = d 4 pi d 4 p2 d 4 P3, d 12 p' = d^Pi d A p' 2 d 4 p' 3 , we have 

X)(^T)*n d 12 p' = J2 ^z'aV d 12 P > (11) 

whenever T and Q are square integrable in terms of the volume element d 12 p'. 
Owing to (10), H and V have the same property [9]. 
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The results of Ref. [7] are as follows: 

a) Provided the three masses are not too much different one from another, equa- 
tions (7) (8) can be inverted in closed analytic form. 

Indeed our model is reliable insofar as, in the no-interaction limit obtained by 
putting the potential equal to zero, one recovers the free motion of independent 
particles. The discussion of this point in Ref. [7] led to require 

\"2 + ^3| < ^H™a, \"2 ~Vz\ < \H m l ( 12 ) 

For instance, if only two masses are equal, say 777,2 = 771,3, their value is allowed to 
deviate from mi by an amount of almost 6%. 

b) Conditions (7) amount to redefine relative energies, in such a way that the 
new relative energies can be eliminated. 

c) The compatibility conditions can be satisfied easily in terms of our new vari- 
ables. Actually, in view of the compatibility requirement, a closed form of the 
interaction is available only in terms of the new variables. 

A typical example of a potential satisfying compatibility and Poincare invariance 
would be of the form 

V = f((z' 2 ) 2 ,(z' 3 ) 2 ,z 2 -z' 3 ,P 2 ) (13) 

since any function of z' B ,y' c and P commutes with y' A ■ P. In this equation all the 
arguments of / are mutually commuting [10]. 

This situation is in favor of using ^ and the new variables, as we shall do hereafter. 
By our transformations, the difference equations (6) become 

4 2 

y' A ■ P * = (-v A - -u B )c 2 a+b (14) 

and the dynamical equation of motion is mapped to 

(3 ^ rnV -P 2 ) * = (D + 6P 2 E + 18V) * (15) 

In order to handle this equation we need to express S in terms of the new variables. 
Lenghty but elementary manipulations reported in [7] show that 

(16) 



(17) 

v (18) 

m + \vl) (19) 
y 3 + \yl) (20) 



s = e + v 2 + 7?e 

where £, 77 are determined by the system 

2.1 . 77 2 

u, v being determined as follows 

P 2 u = y' 2 ■ P + ^y' 3 ■ P - (y 2 ■ 

P 2 v = y' 3 -P + \y' 2 -P- (y 2 ■ 
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2 Three-dimensional reduction 



Now the dependence of on the new relative energies is easily factorized out, 
provided we assume a sharp linear momentum, say 

P a y = k a ^, k 2 = M 2 c 2 (21) 

for some constant timelike vector k. Let ~ denote the projection orthogonal to k. 

^ z • k 

For instance the transverse piece of z with respect to k is z = z ro~k, etc. In the 

k 

rest frame we have y\ = —y\, z\ = — z\, etc. 

We make this convention that, in any operator F depending on the dynamical 

4 2 

variables, the underline indicates that we replace y' A ■ P by {t^a — -^ v b)c 2 and P a 
by k a , hence P 2 by M 2 c 2 . Let us write this symbolically 

F = subs.(y' A -P=(^ A -^ B )c 2 , P a = k a , F) 

It is clear that F reduces to F_ on the " mass- momentum shell" (defined as the 

subspace of solutions to the mass-shell constraints which are eigenstates of total 

linear momentum). 

Note that E depends only on y2,y~3- 

Equation (15) yields the reduced equation 

(3 ™ 2 ~ M 2 )^ = (6(f/i + Vl + m ■ m) + 18Z + 6M 2 c 2 H) $ (22) 

for a reduced wave function ip which depends on three-dimensional arguments only 
(say yi , 2/3 in the momentum representation) . 



3 Equal Masses 

Fortunately, in the case of three equal masses, say m a = m we have this further 
simplification that va = 0, which finally renders u,v of the order of 1/c 2 . More 
precisely (19) (20) entail 

M 2 c 2 u = -{y 2 -m + \yl) (23) 

M 2 c 2 v = -(y 2 -y 3 + \yi) (24) 

hence u and v to be inserted into the reduced version of the system (17)(18). Solving 
for £, 77 we obtain 

I = 2u-v + 0{l/c A ) (25) 

■q = 2v -n + 0(l/c 4 ) (26) 

Inserting this into (16) we get 

E = 3(u 2 + v 2 -u-v) +0(l/c 6 ) (27) 
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correcting a misprint in the higher-order term of equation (100) of Ref. [7]. Since 
the leading term in E is <9(l/c 4 ), we have defined F by setting 

m 4 c 4 s = r = r (0) + ^r (1) + • • • (28) 

where F^, ■ ■ ■ remain finite when c — > oo. In view of (27) it is clear that 

M 4 c 4 H = 3M 4 c 4 (u 2 + v 2 - u ■ v) + 0{l/c 2 ) 
We compute respectively u 2 , v 2 and Hence 

r (0) = \ {(yl? + (y|) 2 + Km ■ m? + 2(g| + vl) (m ■ m) - vlvi) (29) 

expression valid only for three equal masses (this formula was given in Ref. [7] 
without proof). Note that T is a positive operator and would survive in the absence 
of interaction. 

For three equal masses, equation (22) takes on the form 

(9m 2 - M 2 )c 2 i> = 6(y 2 + yf + y 2 ■ y 3 ) 1> + ^V^J + T i> (30) 

Defining 

6A = (M 2 - 9m 2 )c 2 (31) 
and using the rest frame (where va-Vb = —Ya-Yb) we can write 

AV> = (y| + y§ + y 2 • y 3 ) ^ - 3Z V - 1> (32) 

Naturally I\ ) admits an expression identical to (29) in terms of Y3- 
In spite of being three-dimensional, the reduced equation (32) as it stands, is more 
problematic than an ordinary eigenvalue problem. Even if the interaction doesnot 
depend on the total energy (that is: V doesnot depend on P 2 ) the term P 2 E in 
(15), which has no counterpart in two-body systems and yields M 2 c 2 E in (22), 
brings out some energy dependence. It follows that (32), is not a conventional 
eigenvalue equation: through (31) the operator to be diagonalized depends on its 
own eigenvalue. This complication is by no means a drawback special to our model. 
As emphasized in [11] it plagues most relativistic wave equations; the mathematical 
theory of this situation is rather involved, but fortunately this difficulty can be more 
easily handled in a perturbation scheme, provided the unperturbed equation is not 
energy dependent. 

In the rest of this paper we focus on the first relativistic corrections. Therefore we 
solve (32) after expansion in powers of 1/c 2 , taking (31) into account, say M 2 = 
9m 2 + 6A/c 2 . In principle, the exact analytic expression for F is known, and is 
itself a series in 1/c 2 . In fact the knowledge of F^ is sufficient for our purpose. 
Assuming that A remains finite in the nonrelativistic limit, we select these solutions 
that are in some sense "close to" the nonrelativistic Schroedinger equation obtained 
by dropping 1/c 2 in (32). 
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This development is justified insofar as the velocity of light can be considered as 
large with respect to some velocity formed with help of the physical parameters 
defining the system. Practically, the constituent masses and the coupling constants 
involved in the interaction term must be combined as to form a quantity having the 
dimension of speed. In principle one should check that this "characteristic velocity" 
actually has something to do with the average velocities of the constituent particles 
in the slow motion approximation. 

The legitimity and the limitations of this procedure vary according to the analytic 
shape of the interaction term and must be discussed in each specific case. 

3.1 Post-Galilean Approximation 

Let us start expanding in powers of 1/c 2 . Using the rest frame, and assuming that 

i> = V>(o) + ^(i) + • • • , V = Z( ) + ^Z(i) + • • • (33) 
the zeroth order approximation to (32) yields the nonrelativistic limit 

^oV'(o) - (y! + yi + y2 • y3) ^(o) + 3 Z( ) ^(o) = ° (34) 

Setting 

£ ( o) = ^l, U = --V {0) (35) 

' m to y ' 

equation (34) can be re-written as 

E(o) V'(o) = — (y| + y3 + ys • y?) ^(o) + £^(0) (36) 

which is similar to the Schroedinger equation of a nonrelativistic problem with three 
equal masses ( except perhaps for complications resulting from a possible dependence 
of V on P 2 ). Indeed we consider equal masses, thus to = 2mo where mo is the 
reduced mass of either of particles 2, 3, with respect to particle 1. The first operator 
in the r.h.side is nothing but the kinetic energy for a nonrelativistic system of three 
masses to, when the center-of-mass motion has been separated. 

At the first order in 1/c 2 we can, in the last term of equation (32), replace T which 
depends on M 2 , by rV ), which doesnot. In view of (31), in this last term, we can 
also replace M 2 by 9to 2 . Hence 

AV = (yl + y| + y 2 y 3 - sv- J^Lty (37) 

with T(o) bi-quadratic in y. Inasmuch as V is not energy-dependent, the above 
equation still has the structure of a nonrelativistic eigenvalue problem, and can be 
solved by treating the last term as a perturbation. 

More care is needed for most realistic potentials, for which V depends on P 2 , hence 
V_ depends on M 2 c 2 . Fortunately, in several cases, this dependence is of higher 
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order, so that it can be accounted for by addition of an extra perturbation term, as 
follows. Assuming that V_ is as in (33) we have 

i r 

= (yl + yi + y 2 • y 3 )Y> - 3Z (0) ^ - + 3Z d)) ^ ( 38 ) 

Since we do not go beyond first order, let us write A = A( ) + -^A^) . 
For any nondegenerate level A, we have 

A (1) =- <^5+3Z(i)> (39) 
where the expectation value must be calculated in the unperturbed eigenstate ^(o) • 
Binding energy. 

Now we are in a position to calculate, at first post-Galilean order, the binding 
energy of a bound state. This quantity is usually defined through the (linear) mass 
defect [12]. So let us evaluate M — m = M — 3m. Taylor expansion of (31) yields 

(M - 3m)c 2 = - - — ^ + 0(l/c 4 ) (40) 

(M - 3m)c 2 = ^) + 1 ( V) _ % + 0(1/c 4) ( 4i) 
m cr m 6m d 

which yields the first correction to binding energy. 



3.2 Jacobi's coordinates 

Equation (38) amounts to a nonrelativistic problem, formulated in terms of the 
canonically conjugate variables z' A ,y' B . Before we turn to the harmonic interaction 
it is convenient to introduce Jacobi's coordinates that have the virtue of simplifying 
the expression of the kinetic energy. So we perform a linear change from z' A ,y' B to 
R^jIIb, as follows. 

For three equal masses, the Jacobi coordinates R2,R3 associated with q' 2 ,q 3 , are 
defined by the formulas [13] 

R-2 = q' 2 - q' 3 , Rs = ^= (2qi - q' 2 - q' 3 ) (42) 

in other words 

R 2 = -z' 2 + z' 3 , R 3 = (z' 2 + z 3 ) (43) 

Inverting (43) yields 

z' 2 = ^(V3R 3 - R 2 ) z 3 = ^(V3R 3 + R 2 ) (44) 



8 



Since (43) is a linear transformation, it is easy to determine conjugate momenta, 
say n 2 , n 3 , such that [R 2 ,n 2 ] = [R 3 ,n 3 ] = iS and [R 2 ,n 3 ] = [R 3 ,n 2 ] = 0, etc. 
We find 

11 \/3 . 

n 2 = --y 2 + -y 3 , n 3 = — (y 2 + y 3 ) (45) 

Hence inverse formulaes 

y 2 = -n 2 + n 3 , y3 = n 2 + — n 3 (46) 

In equation (36) kinetic energy was expressed in terms of the heliocentric coordi- 
nates. But with help of (46) we can write 

y| + y§ + Y2 ■ y 3 = n 2 + n 2 (47) 

Now (36) may be re- written in terms of the Jacobi coordinates. For the total kinetic 
energy we have 

E^ = r + 1 ( n 2 + n^) (48) 

^ 2m bm m 

In order to compute the first relativistic corrections we need to evaluate also T( ) i n 
terms of II 2 , II 3 . So we must insert (46) into (29). To this end we can write 

|r (0) = A 2 + B 2 + AC 2 + 2(A + B)C - AB (49) 

with this notation 

A = {y 2 )\ B = (y 3 ) 2 , C = y 2 -y 3 

I 3 



A = (n 2 ) 2 , b = (n 3 ) 2 , c = n 2 n 3 



From (46) we get 



2 1 

^4 = n 2 , - — n 2 • n 3 + -n| 

b = n 2 + -^n 2 • n 3 + In 2 
1 

3" 



c = ^n 2 - n 2 



in other words 



Inserting into (29) we get 



ir (0) = (yi + B) 2 + 4e 2 



lr (0) = (nl) 2 + (nl)' 2 + 2n^ni + 4(n 2 ■ n 3 ) 2 (50) 
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4 The Covariant Harmonic Potential 



In order to test the formalism, it is natural to consider first a toy model, namely 
the harmonic oscillator. 

Harmonic interactions are implemented through the potential 

2 2 2 

V = K12 {q'i - q 2 ) + K23 (q 2 - 9 3 ) + K 13 {q[ - 9 3 ) ( 51 ) 
where K a b are positive coupling constants. 

If, for the sake of simplicity, we assume that all these constants are equal we obtain 

this version 

v = ^{q'a-q' b ) = 2k{(z 2 ) 2 + (z> 3 f - z' 2 ■ z' 3 } (52) 

a<b 

where k is a positive constant. The identity 

E (9^) 2 = 2{(z' 2 f + (z 3 ) 2 " z 2 ■ z> 3 } (53) 

reads, after reduction to the rest frame 

E(^a - q^) 2 = 2 \A + z' 2 - z' 2 • z' 3 ] = A + z' 2 + (z' 2 - z' 3 ) 2 (54) 
With help of (44) we have z' 2 — z'3 = — R 2 and 

/2 ^ /Q"D2 o. /OD . T> 1 "D 2 



Note that 



z 2 = - (3R§ - 2 V / 3R 3 • R 2 + R 2 ; 
z'l = (3R| + 2V3R3 • R 2 + R|) 



2[z' 2 +z' 2 -z' 2 -z' 3 ] = ^(Ri+Ri) 



for all choice of units. Hence we obtain 

E(q'a-q' b ) 2 = ^2 + R!) (55) 

which exhibits the Oq invariance of our potential U. 
Finally equation (37) takes on the form 

av = (n 2 + n 2 )^ + \k (R 2 + R 2 )V - (g^W (56) 

For the moment, let us consider the zeroth-order approximation and divide by m. 
We obtain the Schroedinger equation of a non-relativistic three-body oscillator with 
equal masses, written in Jacobi coordinates R 2 ,R3,II 2 ,Il3 (the SUq invariance of 
the nonrelativistic limit would become manifest if we were to choose an appropriate 
unit of lenght). 

In order to make the contact with textbook notations [13], we may define 

K = ^ 
m 
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The nonrelativistic potential is 

t/ = --z = ^-E(q'a-<H) 



3„ K^, , , n2 



m 



At the zeroth order the ground-state wave function is a Gaussian, as well in the 
coordinate as in the momentum representation. We have better to choose the latter, 
where the operator V is multiplicative. Then the unperturbed ground state is 

Y> (0 ) = = const. eX p{-^=(Il2 + n§)} (57) 

In order to check the validity of expanding in powers of 1/c we observe that the 
1 



quantity — \/3K/m has the dimension of a squared velocity. The velocity obtained 
m v 

by taking its square root is a characteristic of the system and should be reasonably 
small with respect to the speed of light. 

Since V doesnot depend on P 2 it follows that V_ doesnot depend on M 2 , thus 

equation (36) is an eigenvalue problem in the conventional sense. 

With notations (35) we have for the n th level of the unperturbed harmonic oscillator, 



E(p) = sJW/m (3 + n), n = 0, 1,2- - oo (58) 

It is convenient to set 

oo = JzK/m = — V2k (59) 
v m 

Indeed, for the ground state we have in particular 

£7(0) = 3w, A (0 ) = 9>/2k = 3mw (60) 



Let us now consider the first post-Galilean contribution; first order perturbation 
theory applies as usual. 

We focus on the ground state; in order to compute the first correction, we need 
to evaluate the expectation value of T( ) in the state V'(o) • At this stage we observe 
that T( ) is a homogeneous function of fourth degree in the six-dimensional vector 
X = (112,113). It follows that, with obvious notations, a being any constant 

< r (0) >= a~ 4 J e~ x2 r (0) PO d 6 X (61) 

provided that a 2 = -(2k)~ 1 ^ 2 , which corresponds to V"(o) = const. e~ Q2x2 / 2 . There- 
fore it is sufficient to carry out the calculation in the case where a = 1, so let us 
provisionally choose the unit of lenght such that k = 2/9. 

It is convenient to note that <\> = 4>2 4>3 introducing the normalized functions 

<p A = ^ exp(-ini) 
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As an operator 11^ doesnot affect 4>b when B / A. Moreover fa, fa are normalized 
to unity, so we have that 

<(ni) 2 >=<fa,(u 2 A ) 2 <p A > 

For instance we obtain 

< (ni) 2 >= ^ (62) 



< (n|) 2 >= l 4 (63) 



and in the same way 

< (T\l\ 2 >= 

4 

Further we have 

< 0203, n 2 , n| fafa >=< fa,u% fa> < fa,nl fa > 

but we compute easily 

3 

< fa, n 2 . fa >=< fa, n| 03 >= - 

thus 

< n 2 n| >= - A (64) 

Finally, if the U J A are the coordinates of the three- vector 11^, we have that 

(n 2 -n 3 ) 2 = (nini + n 2 n 2 + nini) 2 (65) 

For the sixfold integral 

< (n 2 • n 3 ) 2 >= n~ 3 J(u 2 ■ n 3 )V n 2- n 3d 3 n 2 d 3 n3 (66) 

we find 

< (n 2 • n 3 ) 2 >= \ (67) 

Linear combination of all these results yields, according to (50) 

< r (o) > = ^ = 11 + 1/4 
In view of (61) we now revert to an arbitrary unit of lenght and write 

< r (0) >= ^fn (68) 
Apply this result to equation (39) where Y-{\) is supposed to vanish, 



45 k 

A (i) = — 5~2 69 
8 m z 
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It is interesting to evaluate the relative importance of this correction. For this 
purpose consider the quantity 

AA _ 1 A(i) 
A c 2 A( ) 

In view of (59) (60) (69) we finally obtain 

AA 5 \/2k 5 uj 



A 16 m 2 c 2 48 mc 2 
Remind that uo/m is the square of the characteristic velocity. 



(70) 



Now inserting A( ) and A(i) respectively given by (60) and (69), into equation (31) 
we obtain up to 0(1/ c 6 ) 

M 2 _ 54^ 135k 
m 2 m 2 c 2 4m 4 c 4 

or equivalently 

^-9=l 8 ^-f(^ (72) 
m z mc z 8 mc z 

In principle these formulas permit to calculate M at the first post-Galilean order 

when k and m are given. 

But in practice one may be interested in a naive model of baryon. In this case it 
is natural to fix M ( e.g. the proton mass) and adjust m and k in agreement with 

(71) . The most simple possibility is to choose first the ratio M/m within reasonable 
limits discussed below, then extract k from (71) or alternatively extract — » from 

(72) . In this procedure the choice of M/m must allow for a reasonable value of 

UJ 

the characteristic velocity. More precisely, — ~ must be small enough in order to 

mc z 

justify our first-order treatment. It is clear that the more M/m exceeds 3, the more 
our system is rapid. 

Example: 
M 

If _ = 3.03, solving (72) yields 



m 

UJ 



mc 2 



= 0.0100 



so that the critical velocity is about 10 percent of the velocity of light. If M is the 
proton mass (Mc 2 = 920 MeV) we find mc 2 = 303.6 MeV for the constituent quark 
mass. 



5 Conclusion 

Our basic equations involve a unique interaction term and are tailored for allow- 
ing elimination of the redundant variables implied by manifest covariance. In most 
relevant cases, the interaction term looks as if it were made of two-body contribu- 
tions. In fact, the two-body nature of these contributions is not exact, because the 
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transformation from original coordinates to the reducible ones somehow mixes the 
individual variables. For instance the reducible relative variable z' A does not exactly 
match the cluster (1A), and so on. This situation can be interpreted as due to gen- 
uine three-body forces that we have automatically introduced in order to ensure the 
mutual compatibility of the constraints and the possibility of a 3D-reduction. 
In general, this reduction gives rise to a non conventional eigenvalue problem. 

Starting from manifestly covariant basic equations offers several advantages: 
Conceptually we realize that a general description of the system must exist even 
before we assign a sharp value to P. And before we impose a sharp value to P, 
there is no rest frame available yet (P being just an operator) which seems to discard 
a three-dimensional formulation at this stage. 

Our approach yields a Schroedinger-like equation only at the end; the reduced 
wave equation (32) contains a post-Galilean correction which would hardly be de- 
rived from an a priori 3D-theory. 

Another motivation in favor of constraint dynamics is the fact that the contact with 
quantum field theory is more easy in a covariant framework. Actually, this contact 
has been thoroughly established in the two-body case [14] where constraint dynamics 
inspired an improved way of summing Feynman's diagrams (see the "constraint 
diagrams" exhibited by Jallouli and Sazdjian [6]) . Of course, further work is still 
needed in order to determine if the mass-shell constraints of a three-body model 
also can suggest similar simplifications in the three-body case. 

We performed a systematic expansion in powers of 1/c 2 . At least for three parti- 
cles with equal masses, the non-relativistic limit has familiar features and the first 
post-Galilean formulation is tractable: at this order the reduced wave equation is 
similar to a non-relativistic equation modified by an overall perturbation of kine- 
matic origin, supplemented by an additional term which stems from the possible 
energy dependence of the interaction potential. 

Within this framework it is possible to compute for instance the first relativistic 
correction to the binding energy of three given (equal) masses bound by a given 
interaction. Or alternatively, in a simple naive model like the harmonic oscillator, 
one may determine the free parameters (constituent mass and/or coupling constant) 
in order to fit a fixed value of the ground-state mass. 

Although we focused on the first post-Galilean corrections, let us stress that our 
formalism is ready for use at higher orders, with help of equation (32) where T, or 
equivalently S, must be suitably expanded. 

Future work is needed however, for concrete applications: we plan to implement 
spin, consider the case of unequal masses, and improve the contact with other 
approaches [15]. In particular, it may be interesting to re- visit the BS equation in 
terms of the reducible variables employed here. 

Finally, our approach seems to be more specially designed for confined systems. 
When nonconfining forces are present, the occurence of scattering states may rise 
the question of cluster separability which is not addressed here. But even so, our 
picture may provide, if not a complete theory, at least a reasonable effective model 
valid in the sector of bound states. 
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